function [P,T] = generate_PT(?)